Dissipation, noise and vacuum decay in quantum field theory 
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We study the process of vacuum decay in quantum field theory focusing on the stochastic aspects 
of the interaction between long and short-wavelength modes. This interaction results in a diffusive 
behavior of the reduced Wigner function describing the state of the long-wavelength modes, and 
, thereby to a finite activation rate even at zero temperature. This effect can make a substantial 

contribution to the total decay rate. 
In this Letter we shall investigate how stochastic interactions between long and short-wavelength modes affect 
vacuum decay in scalar quantum field theory. Two of the present authors have already considered the relevance of 
stochasticity in the context of the creation from nothing of the Universe Our analysis of that problem led to 

■ the conclusion that the noise- induced transition amplitude was actually larger than the usual quantum estimates [Q . 
However, it remained unclear whether the relevance of stochasticity for the full decay amplitude was a peculiarity 
of gravitationally bound systems, or rather a generic feature of vacuum decay in field theory. The results we shall 
discuss here point quite conclusively in the second direction. In pursuit of clarity, we shall omit most of the technical 

■ details, which shall be reported in separate publications 

ly-j \ As a simple nongravitational example, let us consider a self-interacting scalar field <E> in Minkowski spacetime. The 
f*^; ■ classical action is 

Sren [*] =~J #X (d* ■ 5$ + Af 2 $ 2 - Ig<f> 3 j . (1) 

O ; 

Although we keep h explicit we set c = 1, M has units of lenght' 1 , <I> has units of M\fh and g of M/Vh. For 
ij simplicity, we shall assume that renormalization has already been carried out and that Eq. (|l]) is a good description 
of the relevant dynamics. This means that the parameters M 2 and g may well be renormalization point dependent; 
in any case, any such dependence will be taken as given. This model may be considered as the limiting case of the 
class of models studied by Baacke and Kiselev [|| , when the coupling in the quartic self-interaction is very small. 

We are concerned with situations where the potential displays a local minimum, separated from the absolute 
minimum by a potential barrier. A system of few degrees of freedom, prepared in a false vacuum state within a 
potential well, may decay in essentially two different ways, namely by tunnel effect, that is, going through the barrier 
in a classically forbidden trajectory, or else, by activation, that is, jumping over the barrier (^JtJ. In systems with 
few degrees of freedom, there must be an external agent, typically a thermal source, for activation to be possible. 
Activation results from the system being driven by noise originating in the source. 

In either case, the decay probability follows the Arrhenius law P ~ Ae~ B . In the tunnel effect, B = Se/% where 
h is Planck's constant and Se is the action for the trajectory which goes under the barrier in Euclidean time ||. In 
activation, B = Vg/ksT, where fc# is Boltzmann's constant, T is the temperature, and V s is the height of the free 
energy barrier measured from the false vacuum j|. We can see that activation disappears as T — > 0. 

Our thesis is that in field theories there is a phenomenon similar to activation, even in the absence of an external 
environment and, most remarkably, that this simulated activation contributes to vacuum decay probability even at 
zero temperature. This phenomenon exists because, while vacuum decay concerns mainly the long-wavelength modes 
in the field, these modes evolve in the environment provided by the short-wavelength ones. Due to the time dependence 
of the long-wavelength modes, even if the short-wavelength modes were initially prepared in their vacuum states, these 
will evolve into coherent superpositions of many particle states. The energy to create these particles is provided by the 
long-wavelength modes. On the other hand, it is not possible to predict the exact number of particles to be created. 
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For Bose-Einstein statistics, for example, if N particles are created in the mean, then the dispersion in this number 
is of order y/N (N + 1), and it is never negligible. 

Therefore, we find a dissipative term in the dynamics of the long-wavelength modes, representing the energy transfer 
towards the short- wavelength modes, but also a stochastic element, related to the fluctuations in the energy flux. These 
two terms are related to each other through the fluctuation-dissipation theorems. We must stress that the presence of 
one of them implies the presence of the other as well. Unlike in Kramers' activation, this environment is intrinsic to 
the system. We should point out that because of this same reason, we are not allowed to prescribe the characteristics 
of noise and dissipation independently of the system dynamics. This means that it is not possible in general to assume 
ohmic dissipation or white noise . 

More concretely, tunneling occurs in models where the system may be trapped into a mctastable state, which is 
separated from the basin of attraction of the true vacuum by a potential barrier. There is a saddle point on this 
barrier, representingthe critical bubble, and most of the tunneling dynamics is concerned with motion along the most 
likely escape path [|lTj, which goes through the saddle in the direction of steepest descent. It is possible to identify a 
few degrees of freedom which parametrize the different configurations on this path; the remaining (infinite) degrees of 
freedom describe deviations from the most likely escape route. In the conventional approaches to tunneling, the role 
of these fluctuations is downplayed: they renormalize the action for the few distinguished parameters, and provide a 
prefactor which, after the contributions to the effective potential have been included in the exponential, is of order 
one We intend to focus on the back-reaction of the fluctuations away from the most likely escape path on the 
quantum dynamics of the critical bubble; for this purpose, we shall borrow tools from quantum open systems theory, 
by considering the few distinguished degrees of freedom as a system interacting with the environment provided by the 
transversal fluctuations. 

The technical complexity of the problem increases sharply with the amount of information one wishes to retain 
within the system. For example, one may parametrize the most likely escape path following Ref. JT^] . As shown in 
this reference (and earlier in Ref. [^4|), a time dependent bubble excites the degrees of freedom in the environment, 
thus setting up the kind of dynamical interplay we wish to analyze (see also Ref. |l5|]). To be able to focus on the new 
(dynamical) aspects of the problem, over and above its geometrical aspects, we shall adopt an intentionally simplified 
parametrization of the most likely escape path which retains the essentials of the physics involved. 

Let us return to the scalar field theory above. The potential V [<fi] = ^M 2 c/) 2 — ^gcf> 3 has a stable fixed point at 
4> = and an unstable fixed point at <f) = <f) s = 2g~ 1 M 2 . The former corresponds to zero energy, and the latter to 
E = E s = VM 2 (p 2 /6 in a volume V. For intermediate energies, we may have bound and unbound states. They are 
separated by a potential barrier, which at zero energy extends from <fi = to <p = <p ex it = 30 s /2. To identify the 
relevant modes, we observe that, if we consider fluctuations around the unstable fixed point </> s , then modes with 
wavenumber k > M are stable. The relevant modes, which partake in the tunneling process, have k < M E3j. We 
therefore write the field as <f> = <p + <p, where the first term contains only modes with k < M, and the second term 
contains the short wavelengths; <j) shall be our system. 

In other words, the field represents the average of the full field $ over volumes of order M~ 3 . By construction, <p 
is slowly varying in space; it is technically simplest to handle it as if it were actually spatially homogeneous. We shall 
therefore regard the configurations along the most likely escape path as a sequence of "top hat" field configurations, 
parametrized by a single degree of freedom <fi(t), representing the field amplitude within a domain of size M _1 , 
outside of which the system field vanishes. The center of mass coordinates of the "hat" may be treated as collective 
coordinates in the usual way, and do not affect our results ||l7| . 

If the quantum state of the full field is described by a density matrix p(c/),if,(f)',ip',t), the state of the 4> field is 
described by the reduced density matrix p r (<fi,(j)',t) — J dip p(4>,ip,<fi' ,ip,t) , or equivalently by the reduced Wigner 
function 



f(<t>,P,t) 



To second order in g and leading order in ft, / evolves according to 



{"-/} + nr(r/ + »W/})- y —iA, (3) 



where the curly brackets are Poisson brackets, H s ~ </> 2 /2 + V((f>), T — J dt' H (t — t') <fi (t 1 ) , and N — 
J dt'N (t — t') 4> (f) . The kernels H and N represent the effects of dissipation and noise, respectively. They come 
from the quadratic part of the Feynman- Vernon influence action |T§| ]. Computing the influence functional requires 
handling formally infinite quantities (and in our case, also a linear term in 0); regularization and renormalization 
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leave a finite residuum, which are the one-loop correction to the effective potential and a finite wave function renor- 
malization. As we have already remarked, we assume that these corrections are already included in Eq. (|l|). This is 
sensible because, although they may be quantitatively important, they do not affect the nature of the problem |Q. 

Our approach to Eq. (||) will be the following. It is clear that if only the first term of the right-hand side is kept, 
the equation reduces to the classical transport equation and there is no tunneling. Retaining the first and last term on 
the right-hand side is equivalent to writing a Schrodinger equation for the wave function of the homogeneous mode, 
as if it were a closed system. Since this is a one-dimensional problem, the tunneling rate may be computed either by 
the instanton or the WKB method, which are known to be equivalent in this case. We wish to know if the middle 
(back-reaction) term makes a substantial contribution to the total rate. With this strategy in mind we shall discard 
the third contribution to the right-hand side in Eq. (^) , assuming implicitly that the back-reaction term is dominant 
(see below). It then reduces to Kramers' equation which may be seen |0j to describe the evolution of an ensemble of 
points evolving according to the Langevin equation: d(f>/dt — p(t), 

A j- t (t) = -V [d> (t)] +jdt'H(t- t>) 4, (f ) + & (4) 

with initial conditions (<fii,Pi) weighted according to the initial Wigner function, and Gaussian noise characterized 
by (£(£)£(£')) = frN{t — t'). Although this representation of the dynamics has an important heuristic value, only 
the reduced Wigner function / has a direct physical meaning, and it will not allow an interpretation as a classical 
distribution function in general (it will not be generally positive definite). 

We are interested in the weak dissipation limit, as discussed by Kramers ||, when the relaxation time is long 
compared with the classical period of motion. On the other hand, the memory time in the integrals in Eq. ([5]) is 
determined by the frequencies in the environment, which are large with respect to the dominant frequencies in the 
system. Thus we are allowed to (and, in a formal expansion in powers of H, we must) use solutions to the classical 
equations of motion within the memory terms; these solutions may be written down explicitly in terms of elliptic 
functions. This is a less drastic approximation than the Markovian one discussed in Ref. |l9) . We also neglect transient 
terms (or in other words, we assume t 3> M^ 1 ); this means that we can take the lower limit of the time integrals in 
Eqs. M) as t — — 00. In the weak dissipation limit the reduced Wigner function depends only on the action variable 
,/ = y~ § d<fi p, and, averaging over angles, Kramers' equation reduces to a one-dimensional Fokker-Planck equation 
df jdt + d<t/dJ — 0, where $ = —<dfl~ 1 df/dJ — A/ is the flux, and f2 = £1 (J) is the frequency of the corresponding 
classical motion [^Cj. The point of this analysis is that it is possible to derive explicit expressions for the coefficients 
O and A in the Fokker-Planck equation Near the value J s of the action variable at the separatrix (that is, the 

limiting trajectory which connects to the unstable equilibrium point) these have finite values, while as J — ► they 
go to zero as E 3 , where E = H s (J). In the thermal activation problem one finds an identical equation, but the 
coefficients decay linearly on E || . 

The weakness of noise and dissipation in our (vacuum decay) problem reflects the origin of these effects in particle 
creation. Since particles are created in pairs, there is a threshold for particle creation at frequency uj ~ 2M. At low 
energy, classical motion is mainly harmonic with frequency M for small oscillations around the mctastable minimum, 
hence particle creation is weak. It never actually vanishes, though, because at any finite energy there is a small 
deviation from harmonic motion. The amplitude of the component with frequency nf2 decays as E n as E — *■ 0, which 
is enough to trigger particle creation [ fH| ]. 

The Fokker-Planck equation describes an initial value problem subject to nontrivial boundary conditions at J = 
and J = J s . These are vanishing flux $ = at J = 0, and vanishing probability / = at J = J s . The linear 
operator L which is defined by Lf = d^/dJ is sclf-adjoint with respect to an adequate inner product ||, and the 
equation may be solved by an expansion in normal modes in the usual way. A general solution is reconstructed as 
a superposition of modes f r decaying as exp (— rt). For a given r, f r oscillates as J — > 0, and the modes must be 
subject to a continuum normalization, as in the usual treatment of the WKB wave function in quantum mechanics 
p2| . The result is that, given any smooth initial condition with mean energies of the order of the false vacuum 
energy HM/2, the persistency amplitude P (t) = 2ir J dJ f (J,t) decays exponentially with a constant A for \t > 1, 
turning to 1/t for longer times (this crossover is also observed in the usual tunneling amplitude p3| ). The constant 

is H A » Aexp {— J dE A (_E)} ~ Aexp |— a-^p- j, where A is of order 1, and a ~ 0.2. By contrast, the tunneling 
amplitude, in the corresponding approximation of only considering the homogeneous mode, yields a similar formula, 
but with a - 4.8 g. We can see that, in this case, the zero temperature activation rate is higher than the tunneling 
amplitude by an order of magnitude in the exponent. 

We should point out that the fact that in our example the activation amplitude is actually larger than the tunneling 
amplitude is model dependent. Roughly speaking, low and broad barriers favor activation, while high and narrow 



3 



barriers favor tunneling. We must also stress that these results must be considered as preliminary, pending a more 
satisfactory parametrization of the system, and therefore a more realistic modelization of the system-bath interaction. 
It is safe to conclude, however, that activation should not be discarded a priori, but rather counted on as a potentially 
significant contribution to the overall decay amplitude. 

Note that we reach a different conclusion to that of Ref. [p4j, where dissipation suppresses tunneling. However a 
direct comparison between the two analysis cannot be done as the two models differ essentially in the coupling with 
the environment degrees of freedom, which is quadratic in our case. Our analysis is closer to that of Ref. |25| ], and it 
is certainly compatible with their conclusions; see also Ref. pfj. 

Still, we must stress that we should not expect a similar behavior in systems with few degrees of freedom. The fact 
that in our problem the environment actually contained a large enough number of degrees of freedom as to represent 
a continuum for all practical purposes is essential to provide a suitable driving force. If some frequency intervals 
were lacking, then there would arise islands of stability where no resonance is strong enough to move the system 
forward. These islands would act as absolute barriers to noise-induced decay, or at least would depress the noise 
induced amplitude much below the tunneling estimates. 

In conclusion, we have shown that vacuum decay in field theory is qualitatively different from the same process 
in systems with few degrees of freedom, because the former are intrinsically open systems. Interaction between 
long and short-wavelength modes induce a stochastic dynamics for the former and results in activation even at zero 
temperature. 
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